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We considered the phantom cosmology with a lagrangian L — — + -qg^'-'^^'^ ^0. y] — 

which is original from the nonlinear Born-Infeld type scalar field with the lagrangian 
L = —[1 — -v/l — rigf^^cl)^ y] — u{4>). This cosmological model can explain the accelerated 
^ ■ expansion of the universe with the equation of state parameter w < —1. We get a sufficient 

I condition for a arbitrary potential to admit a late time attractor solution: the value of 

OO , potential u{Xc) at the critical point (Xc,0) should be maximum and large than zero. We 

^\ . study a specific potential with the form of u{(f)) — Vo(l + *o via phase plane analysis 

o ■ .... 

' and compute the cosmological evolution by numerical analysis in detail. The result shows 

' that the phantom field survive till today (to account for the observed late time accelerated 

expansion) without interfering with the nucleosynthesis of the standard model(the density 
parameter fl^ ~ 10~^^ at the equipartition epoch), and also avoid the future collapse of the 
i-C ' universe. 
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1 Introduction 



A recent great finding in modern cosmology is that the universe is undergoing a accel- 
erating expansion at present epoch and have a decelerated expansion phase in the recent 
past. This dynamics property is not only inferred in high redshift surveys of type la super- 
novae[l], but also independently implied from seven cosmic microwave background experi- 
ments (including the latest WMAP results) [2]. People have proposed many candidates for 
dark energy to fit the current observations, several important candidates may be the cos- 
mological constant (or vacuum density), quintessence[3](a time varying scalar field evolving 
in a specific potential) and the "k-essence" (see literatures[4,5] and References). The major 
difference among these models is that they predict different equations of state parameter of 
dark energy and different accelerating style of present universe. Very recently, joint analysis 
of CMB-|-SN-Ia-|-HST+LSS data hold that the equation of state parameter can lie in the 
rang of —1.32 < ti; < —0.82 with a best-fit value of it; ~ —1.04 at the 2 — cr confidence levels, 
slightly preferring "phantom" models [6]. Though phantom models [7] poses a challenge to the 
wildly accepted energy condition [8, 9] and may lead to rapid vacuum decay [10], it is still very 
meaningful to study these models in the sense that it fits the observation very well, and also 
is phenomenologically interesting. 

Nonlinear Born-Infeld theory has been considered widely in string theory and cosmology. 
In 1934[11], Born and Infeld put forward a theory of nonlinear electromagnetic field to resolve 
the singularity in classical electromagnetic dynamics. The lagrangian density is 



In order to describe the process of meson multiple production connected with strong field 
regime[12], Heisenberg proposed the following nonlinear scalar field lagrangian firstly: 



This lagrangian density(2) possesses some interesting characteristics: (i)it is exceptional in 
the sense that shock waves do not develop under smooth or continuous initial conditions[13], 
(ii)because nonlinearities have been introduced, nonsingular scalar field solutions can be gen- 




(1) 




(2) 



erated, (iii) if g'^'^f/)^ ^ << ^, by Taylor expansion, Eq.(2) approximates to the lagrangian 
of linear scalar field. 
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the linear theory is recovered. H.P.de Oliveira has investigated quahtatively the static and 
spherically symmetric solutions of this nonlinear scalar field [14]. One of our authors has 
investigated the quantum cosmology based on the same nonlinear B-I type scalar field[15]. 
Especially, if the potential u{(f)) equals to |, we find that pp = this describes a Chaplygin 
gas. The Chaplygin gas has raised a renewed interest recently because of its many remarkable 
and intriguing unique features[16]. 

We have proposed a dark energy model based on the lagrangian Eq(2) in literature [8], 
and found the universe with this B-I type scalar field with a potential can undergo a phase 
of accelerating expansion, the corresponding equation of state parameter lies in the range of 
— l<w<—^, for the nonlinear B-I lagrangian with a negative kinetic energy term, when 
u{(p) is constant, the corresponding equation of state parameter can lie in the range w < —1. 
In this paper, we consider the phantom cosmology of this nonlinear B-I lagrangian, i.e, with 
a negative kinetic energy term, L = ^ [1 — ^yl + rjg^^'^cj)^ ^(f)^ u ] — u{(f)), where u{(j)) is taken 



the form of n(<^) = Vb(l + ^)e '*o , investigate the global structure of the dynamical system 
via phase plane analysis, and calculate the cosmological evolution by numerical analysis. 
The paper is organized as follows: In section 2, we start the phantom model with the B-I 
lagrangian, in section 3, the dynamical evolution of the phantom field without the presence 
of matter and radiation is considered, we discuss the evolution of the phantom field at the 
presence of matter and radiation in section 4, section 5 is summary. 

2 The Phantom Model with B-I Lagrangian 

We consider the kinetic energy term is negative, namely the phantom cosmology. For the 
spatially homogeneous phantom scalar field, we have the following lagrangian: 



where u{(f)) is the positive potential, r/ is a constant. First we consider the simple case 
where the phantom field (j) is dominant only. In the spatially fiat Robertson- Walker metric 
ds^ = dt^ — a^{t){dx^ + d'lp' -|- dz^), the Einstein equation G^v = ^T^i/, can be written as 




(4) 




(5) 
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For a spatially homogenous phantom field (p, we have the equation of motion 

4> + 3H^{l + ri^'^)-u{(t)){l+ri^'^)i =0 (6) 

where the overdot represents the differentiation with respect to t and the prime denotes the 
differentiation with respect to 0, the constant ^ = SttG, where G is a Newtonian gravitation 
constant. The Energy-moment tensor is 

r'^ = -£^£l^tl^=-S^'L (7) 



From Eq.(7), we have 



r/y 1 + rj^i^ ' 



1 V 1 + rid)^ 
T] ri 

From Eqs.(8) and (9), we obtain 

P^+P<t> = — , (10) 

It is clear that the equation of state parameter Wtf, < — 1 is completely confirmed by Eq.(lO) 
which agrees with the recent analysis of observation data. We also can get 

i2 



p</, + 3p^ = - - Vl + r?(^2_ _2u (11) 

"I V 1 + 

It is obvious that + 3p^ < 0. Eq.(ll)shows that the universe is undergoing a phase of 
accelerating expansion. The model of phantom B-I scalar field without potential u{(j)) is hard 
to understand. We can always find = — ■j^=^= < in this case, it is unreasonable 

apparently. However, in the model of phantom B-I scalar field with potential u{(j)), when 

<4>) >l- I ^ . (12) 

Pff) is always greater than zero. In phantom B-I scalar model with a potential t((0), we also 
find the strong and weak energy condition always failed from Eqs.(10)and (11). The equation 
of state is 

= — = -1 , (13) 

P't> 1 + [r]u{(t)) - I + r]^'^ 

It is always less than —1 if Eq(12) is satisfied. 



We also can write the speed of sound in these phantom model: 



X 



La 



X 



P, X ^ X + 2XL^ ^ X X 

where P = L^{(j), X) and p = 2XL^^x — L(j,{(p, X) with X = 5(9^^)^. In our model X 
so the speed of sound of the phantom with BI lagrangian is 

2 Px 



P, X 



= 1 + ry^^ 



(14) 



(15) 



It is clear from Eq.(15) that the Cg is always greater than 1 unless = 0, this means that the 
perturbation of the background field can travel faster than light as measured in the preferred 
frame where the background field is homogeneous. For a time dependent background field, 
this is not a Lorentz invariant state. However, it does not violate causality because the 
underlying theory is manifestly Lorentz invariant and it is not possible to transmit information 
faster than light along arbitrary spacelike directions or create closed timelike curves [17]. In 
Ref[18], they also shows that a fluid with | li; |> 1 does not constradict causality if w is not 
constant. In our model, w is clear a time-dependent parameter. 



3 Dynamical Evolution Of The Phantom Field 

As we consider the phantom field becomes dominant ,we can neglect the nonrelativitic 
and relativistic components (matter and radiation) in the universe, then from Eq. (5,8,6), we 
have 

(^ + (^(l + #2) 



- + u((/))]-m((/))(1 + #^)2 =0 



(16) 



to study an numerical computation, it is convenient to introduce two independent variables 



X = 
Y = 



(17) 



then Eq.(16) can be written 



dX -V 
dt ~ ^ 

^ = u{X){l + r,Y^)i - ¥{1 + r,Y^)j3^[-j^ - i + u{X)] 



(18) 



we can obtain this system's critical point from 



then its critical point is (Xc,0), where the critical value is determined by u'{Xc) = 0. 
Linearizing Eq.(18) around the critical point, we have 

dX _ Y 

(20) 

^ = u"{x,){x - X,) - /(sRx^Jy 

the types of the critical point are determined by the eigenequation of system 

+ aA + /3 = (21) 
where a = ^2,iu{X^),f3 = -u"{X,), the two eigenvalues arc Ai = -V3^"(^^)+\/3g"(^c)W(Xe) ^ 
As = - V3gn(x,)+4n"(x,) _ ^ positive potentials, if u"{X^) < 0, then the criti- 
cal point (Xc,0) is a stable node, which implies that the dynamical system admits attrac- 
tor solutions. We can also conclude that if a potential possesses the general properties: 
u{Xc) > 0,u'{Xc) = and u"{Xc) < 0, then our phantom model with this potential must 
have a attractor solution and will predict a late time de-sitter like behavior (f/;<^ = — 1). Some 
authors have studied the cosmological dynamics of phantom scalar field with special potentials 
which possessed above properties, such as u{(f)) = Voe~^^^'^^[19], u{4>) = Vo[cosh(-^]~-^[20], 
u{4>) = VQ[l + a(t)''/Mp\-\ 

Next we specify the potential with a special form. We choose a widely studied potential[21] 

as 

u{4>) = Vb(l + ^)e^"i^ (22) 

where Fq > 0. It is easy to find that the critical Xc = and u"{Xc) = — ^ < in such a 
potential. Therefore this model has an attractor solution which corresponds to its attractor 
regime, the equation of state w < —1. Substitute Eq.(22) into Eq.(18), we obtain 

dX 

dt ~ ^ 



^ = -|^e(-^)(l + vY')l - Y(l + vY')^3a^^^^ - \ + Vo{l + f^-o^] 

(23) 

We will solve this equations system via the numerical approach, to do this, we rescale the 
quantities as x = s = i4>o'r])~^t, y = ^J^Y ^ Then Eq.(22) becomes 

dx — ^ 

dS " IC\A\ 

^ = -^xe(-)(l + 2/2)1 - <^oy(i + y2y3^[_i_ _ 1 + ^(1 + a;)e(-)] 

where 7 = Vorj are parameters, we take ^ = 1 for convenience. The numerical results with 
different initial condition are plotted in Figs.l — 3 and the parameters (l)o = \/0T, 7 = 3. 



Figl.This plot shows the evolution of the scalar field in difference initial condition, 
solid line is for 0j„ = O.70O, dotted line is for 0j„ = 00, dashed line and dot- 
dashed line for0j„ = 1.3<po, ^■'7(l>o respectively, they are all plotted for a fixed value 
of Vin = 0.1 . 




Fig2. The evolution of w with respect to s, the initial condition is the same as figl. 













/ / / /' \V\\ 


1 1 

] ) 


] j 


-0.5- 






/ / 


y -1- 








-1.5- 






,/'' 


-2- 








-2.5- 









0.5 1 1.5 



Fig3. The attractor property of the system in the phase plane, the initial condition 
is the same as figl. 

As we know, when r/ ^ 0, our model come back quintessence case. In order to see the 
nonlinear effect, we plot the quintessence model with our model in fig4 and fig5. 



X 




1 2 3 4 5 6 



Fig4. the evolution of scalar field with respect to s, solid line is quintessence model, 
dotted line and dashed line are nonlinear scalar field,dotted line is for r/ = 2/3,dashed 
line is for = 1/3. 
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Fig5.the evolution of w with respect to s, solid line is quintessence model, dotted 
line and dashed line are nonlinear scalar field, dotted line is for 77 = 2/3, dashed line 
is for ?7 = 1/3. 

From figl,fig2 and fig3, we can easily find the system admit a attractor solution, the 
equation of state parameter w start from nearly -1, then quickly evolve to the regime of 
smaller than -1, turn back to execute the damped oscillation, and reach to -1 eventually for 
ever. Due to the unusual physical properties in phantom model, the phantom field, released at 
a distance from the origin with a small kinetic energy, moves towards the top of the potential 
and crosses over to the other side and then turns back to execute the damped oscillation 
around the critical point, after a certain time the motion ceases and the phantom field settles 
on the top of the potential permanently to mimic the de Sitter-like behavior (w^ = —1). Fig4 
and fig5 indeed shows that the nonlinear scalar field will come back to quintessence case 
when r] decreases to zero, the nonlinear effect does not affect the global attractor behavior 
but change the evolution of w and scalar field in details. 

4 Evolution of the phantom at the presence of matter and 
radiation 

In the previous section, we studied the case when phantom dominates over all other energy 
density. Now we consider the more general case that the effect of matter and radiation can 
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not be neglected, thus the Eq.(5)becomes 

= |(P0 + m + Pr) (25) 

where pM and pr denote the energy density of nonrelativistic and relativistic components(namely, 
matter and radiation), respectively. We know that 

Pmo^ = pMiO-i = constant 
Pror = priOl = constant 

So we can rewrite Eq.(25) as 

ij2 ^ ^2[Pi + j^^.(^)3 + j^^^(^)4] ^27) 
Pci ^ Oj 

where Hf = pd is the critical 6n6rgy d.6nsity of th.6 univ6rs6 Eit th.6 initicil ti , /^m^ , /^t-^ cLr6 
the matter and the radiation enerffy densities at t-i ^ J^^^ ^ ^t^? 

are the cosmic density parameters 
for matter and radiation at t^, ai is the initial scale factor at ti, we will specify = 1 for 
convenience. We can finally give the form of Eq. (27) as 

^2 = H^[P± + n^.a-^ + n^.a-^] (28) 

Pci 

we introduce the new variables x = = = ^i^^.^ we can express Eq.(16) as 



(29) 



5l = -3|/(i + y^) + V^^«^V' 



where ip{x, y, N) = [-!- ( — ^ . - ^ + «) + ^Mie + J^r-iC ^■'^] 2 . Substitute the potential 



Eq.(22) into Eq.(29), we get 

dx ^ a ^/ 

[ ^ = -3y(l + y^)-^x(l + y^)ie-V 

where V'(x,y, A^) = i^i-^^^^ ' Wo + + ^)^^""^) + ^Mie'^^ + J^Mie'^^J-i The 
numerical results are plotted in Figs. 6 — 10. We specify our starting point as the equipartition 
epoch, at which D,Mi = ^ri = 0.5, We plot the results by choosing the parameter ^ = 10~^. 
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Fig6. The evolution of the scalar field with respect to N in the presence of mat- 
ter and radiation, they are plotted in different initial conditions, solid line is for 
(f)in = 10.6^0 jdotted line ,dashed line and dot-dashed line for 0j„ = 8^0)600)400 
respectively, yi„ is set zero all the time. 
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Fig7. The evolution of w with respect to N in the presence of matter and radiation, 
we chose the same initial condition as fig6. 
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Fig8. The attractor property of the scalar field in the presence of matter and 
radiation, the initial condition is the same as in fig6. 
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Fig9. The evolution of density parameter $7 with respect to N , solid line is for scalar 
field ,dotted line for matter and dashed line for radiation, where the initial value of 

(pin = W.GcpcUin = 0. 
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FiglO. The evolution of decelerating factor q with respect to N, the initial condition 

is the same as fig9. 

From these results(fig6,7,8) we can conclude that there is still a late time attractor solution 
even in the presence of matter and radiation. The phantom field energy density will domi- 
nant over the matter and radiation composition(J70 — ^ l)with its state parameter — — 1, 
mimicking the cosmological constant. Such evolution behaviors will avoid the cosmic dooms- 
day [9]. FiglO quantitatively shows the universe is speeding up at present but undergoes a 
decelerated regime at the recent past. Since Tgq c± h.QA{p,Qh?)eV 2± 2.843 x 10^ii'(where we 
set h = 0.71), To ~ 2.7K, ai = Ueq = 1, the scale factor at the present epoch oq would nearly 
be 1.053 X 10^, then we know at present epoch N = Ina^ ~ 9.262, we calculate the current 
value of r20(, ~ 0.705, decelerating factor qq ~ —0.603 and the transition shift Zt — 0.672 in 
our model, this is consistent with current observation data. 

5 Summary 

We have shown that it is possible to construct a viable cosmological model based on 
nonlinear Born-Infeld scalar field. We get the condition to admit a late-time attractor so- 
lution: u{Xc) > 0,u'{Xc) = and u"{Xc) < at critical point, this condition is also true 
for quintessence model. However, we should emphasized that this condition is a sufficient 
condition. The attractive behavior of the solution may resolve the "coincidence problem"' in 
dark energy models in principle. We investigate the global structure of the dynamical system 
via phase plane analysis and calculate the cosmological evolution by numerical analysis with 
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and without the presence of radiation and matter, the result shows the universe wiU evolve 
to a de Sitter like attractor regime in the future and the phantom field energy density will 
dominant whole universe and behave as a cosmological constant. An obvious excellent char- 
acter of our model is that the phantom field survive till today (to account for the observed 
late time accelerated expansion) without interfering with the nucleosynthesis of the standard 
model(the density parameter $7^ ~ 10^ at the equipartition epoch), and also avoid the 
future collapse of the universe, it was indicated by Carroll, Hoffman, and Trodden [19] 

Theoretical cosmology with phantom models has become an active area of theoretical 
research. However, recently, some authors showed in [22] that quantum effects (without a 
negative energy kinetic term) can also render w < — 1 on cosmological scales. Since the idea 
of phantom cosmology is pretty new, there are many questions remaining open yet. 
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